I1HAC0, Inc. 735 W. Clinton Street Ithaca, New tor-: 14 850 

-V " • 3 L- .* y'_ jj. 

February 22, 1972 
Report No. 90345 
File No. 10-2460 
Approval 

ABSTRACT 

A Magnetic Control System for Attitude Acquisition 



This report describes a spacecraft magnetic attitude 
acquisition system that is capable of automatically despinning 
a satellite from arbitrarily high rates arcund any axis, and 
provides terminal orientation that makes capture by conventional 
fine c:- l'rcl attitude control systems routine. The system con- 
sists ci a 3- axis magnetometer, a set of 3 orthogonal magnets, 

•ju.d appropriate control logic. No earth sensor is required. 

Acquisition is treated in two phases. During the despin phase 
we are concerned v.ith removing the tumbling motion of the satellite 
In this phase the performance of the system is unaffected by the 
presence or absence of a momentum wheel. Phase 2, that of 
orienting the spacecraft to the desired attitude, requires that a 
m oment um bias wheel he present. In the terminal orientation, the 
&cis of' the mOliien c Unfwh eel is Substantially perpendicular to the 
orbit plane (roll and yaw errors near zero) and the pitch rate of 
the satellite is at twice orbit rate. 

This report describes the analysis and simulation that has 
been done in evaluating the performance of this system. A well- 
configured system will result xn despin times of the order of $ 
orbits per RPM for spacecraft in low earth orbits. Following 
despin, terminal orientation is achieved after another one to 
three orbits, depending on the capture range of the associated 
fine control system. 

While this report does not describe the physical hardware, 
the system can be implemented inexpensively v.ith weight less than 
5 lbs. and power of about 3 watts. 
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INTRODUCTION 

In this report we will be concerned with the description and 
evaluation of a magnetic acquisition/attitude control system. The 
work described here has been performed during the last several 
months, principally by the author, with the assistance and guidance 
of Messrs. David Sonnabend, and Robert Fowler. This work was done 
under contract #NAS5~21649, at \he behest of the Systems Analysis 
Branch, GSFC, Greenbelt, Maryland. The original concept of this 
system was suggested by Seymor Kant, Head, Systems Analysis Branch, 
GSFC. 


The attitude control system described here is designed to 
"despin" a satellite and then orient it to any preferred Earth 
referenced attitude. This is done thru interaction with the 
Earth's magnetic field. An interesting feature of the system is 
that it requires no Earth sensor. 

In this report we will discuss the control law for despinning 
the sate lite and the mechanism thru which some preferred attitude 
is attained. Then we will discuss stability considerations as they 
apply to this system. Next, system performance is evaluated in the 
light of both analytical results and simulation results. We will 
be concerned with both despin and attitude acquisition performance. 
Finally, some recommendations for further study are included in 
case any parties might wish to continue the work discussed here. 

Appendix B displays a listing of the FORTRAN program used in 
the digital simulation. Appendix C reviews some of the dynamic 
stability theory for dual spin satellites, for convenient reference, 
but does so in a different fashion than is usually employed. The 
results obtained are, of course, the same as usual. In Appendix D 
we analyze the equilibrium points of our system; the stability of 
these points is briefly treated in Appendix E, using some of the 
results of Appendix C. 

The author wishes to make one further remark before continuing. 
The (differential) equations describing the behavior of the system 
discussed herein are nonlinear and otherwise intractable. This 
being the case, it is not possible to obtain explicit analytical 
solutions and recourse to simulation must be made. Much of the re- 
sults presented here were obtained via a digital simulation model. 

A program, previously written by the author and noted in the refer- 
ences (1) , was used to this end. A listing of this program, and an 
index of computer runs, may be found in the appendices. A brief 
description of the simulated satellite is contained in the following 
section. 


1 
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SIMULATED SATELLITE 

Most of the simulation runs dealt with a satellite with the 
following characteristics: (1) moments of inertia of 47.5Kg-m 2 

(35 slug- ft 2 ) , 67.9Kg-m2 (50 slug- ft 2 ) , and 33.9Kg-m 2 (25 slug-ft-) 
about, respectively, its yaw, roll and pitch axes; (2) a pitch 
momentum bias of 0.944 Nt-m-sec (0.7 lbf-ft-sec) ; and (3) switch- 
able control magnets of 10 4 pole-cm strength. There are three such 
magnets, one each along the yaw, roll and pitch axes. We again 
note that these are principal axes for the satellite. This satel- 
lite is in a circular low altitude orbit with a period of approxi- 
mately 100 min. We note however, that these parameters of orbit 
and vehicle are of no special significance, but were chosen merely 
to implement the digital simulation. 

DESCRIPTION OF PROBLEM 

In general, a satellite injected into orbit will be observed 
to be tumbling about with a certain residual angular velocity. This 
is so in spite of efforts made to release it "gently”. Various de- 
vices have been utilized to dissipate this residual motion. In our 
case, the satellite will be launched by a spin stabilized rocket and 
hence will be inserted into orbit with an angular velocity of roughly 
100 RPM about some (known) axis. A simple mechanism (a yo-yo) will 
then reduce that spin to the order of 1-10 RPM. However, this rate 
is, in general, still so high as to prevent the basic attitude con- 
trol system from locking onto some target (the Earth, sun, stars, 
etc.). A further spin rate reduction is required, preferably to say 
0.1 R*V'i or less. To this end, such active devices as rate gyros, 
accelerometers, and magnetometers coupled to gas jet systems, and 
passive devices such as eddy current and magnetic hysteresis rods 
have been employed. 

Reverting to general terms, then, there are two distinct 
phases of the attitude acquisition procedure. The first is despin, 
of which we have just spoken, and the second is an attitude orien- 
tation phase which facilitates capture by the basic control system. 
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CONTROL THEORY 

As noted above, the attitude acquisition system performs two 
principal functions. The first is despinning the tumbling satel- 
lite, and the second consists of providing some preferred orienta- 
tion. Consider the former. 

The rotational kinetic energy of the satellite, not including 
that of any reaction wheels, is 

2 

T - 1 l (3) 

2 i 

We would like to diminish this quantity. 

d (T) = T = = T * w (4) 

dt i 

_ where t is the external torque acting on the satellite, and 
t • w is the rate at which work is done on the satellite. 

For a dipole M in a magnetic field B, the torque exerted by the 
field on the dipole is 


T = M x B 

Combining (4) & (5) 


(5) 


T = M x B 


to = B X W 


M 


BODY — 

Now consider the quantity d_ (B) 

dt 

d SPACE (B) = d B0DY (B) + w X B 


( 6 ' 


(7) 


dt 


dt 

.SPACE 


(B) =0 . That is, we say that this term is negligible 


Suppose 

dt. 

in comparison with the other two terms. . Then (7) becomes 

d B0DY (B) = B x w (8) 



in (6) we obtain 

M 3 ■ /\ ■ ■ 


/ J 

VAV 




(9) 


’TKiT'clesired control scheme is now clear. We must simply measure 
B along any axis in the satellite, and then change_the polarity of a 
magnet lying along that axis, keeping the sense of M opposite that of 

B. This insures T < 0 and decreasing T. The limitations of such a 
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scheme are also clear. Equation (9) is valid only when (8) is 
valid, that is, if d SPACE (B) is truly negligible in comparison 

_dt 

with w x B. Now, B varies in inertia space_as sin2w 0 t, where 
u>o is the orbital frequency. Thus, d sp ^CE (B) is equivalent to 

_dt 

a body rate of 2u 0 r and as long as w is, say, an order of magni- 
tude greater than 2u) 0 (about 0.002 rad/sec for low altitude 
orbits), we may safely neglect d SPA CE (B) . 

dt 

For reasons which will become clear later, the chosen con- 
trol configuration involves three magnets.- These lie along the 
principal axej of the vehicle and are controlled by magnetometers 
(which measure the field component) along these same axes. The 
controller operates in a flip-flop manner, switching the magnets' 
polarity to keep the M-^ opposite in sense to fei, as measured by 
magnetometers on board the satellite. 
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TERMINAL ORIENTATION 

We now consider the second phase of attitude acquisition, 
orientation. We wish to attain some particular attitude. As long 
as this preferred orientation is with respect to the Earth, any 
desired attitude may be attained. In this section we will confine 
our discussion to a satellite having only a pitch momentum bias. 
Other momentum wheels, if any, are not yet activated at this point. 
The arguments advanced below are not as rigorous as might be de- 
sired but it is felt that their physical appeal more than compen- 
sates for that. i 

Consider the satellite to be orbiting the Earth at angular 
rate u 0 / and assume also that it has been completely despun. Then 
it must be that the satellite interprets this situation as though 
it were spinning backwards at angular velocity -2to 0 e n , where e n is 
a unit vector normal to the plane of orbit (and whose sense is 
determined by the direction of orbit) . This is because the satel- 
lite senses the Earth's magnetic field passing by it twice per 
orbit. Note that the trajectory yaw and roll axis components of 
the Earth's field vary as sin 26, where 0 is the angle into orbit 
from some reference point. Since the system "senses" an angular 
velocity whose direction is -e n , it responds with a corrective 
torque along the direction +e n * The net result of such a torque 
must be to gradually align the system momentum bias (which properly 
lies along the vehicle pitch axis) with the orbit normal. In this 
way the proper yaw- roll attitude is attained. It is the fact that 
the satellite is not spinning at 2wo that results in the net torque 
along e n which in turn aligns the momentum bias with the orbit 
normal. But what if it had been spinning at 2wo originally, or in 
some way reached this state. It is apparent that this could occur 
in only two ways. That is, the system momentum must be either 
parallel or antiparallel to e n . If this is not so, then a component 
oi tnis momentum lies in the orbit plane and a body fixed torque, 
also lying in the orbit plane and of the correct magnitude would be 
required to force the body to rotate about en at 2coo* Where is this 
torque to come from? The author does not believe in its existence 
and hence concludes that the system momentum bias must lie parallel 
or antiparallel to e n if it is to rotate at 2u) 0 . Both of these 
situations correspond to equilibrium points. We simply note here 
that the antiparallel situation is not a stable condition. 

Our conclusions are then that the satellite will eventually 
align its momentum bias axis (in our case, the pitch axis) with the 
orbit normal and thus provide yaw-roll attitude acquisition. If 
the reader is not completely convinced of this by the racher un w 
rigorous arguments advanced here, we reassure him by noting that 
considerable simulation experience is completely in accord with the 
above analysis. A further discussion concerning the amount of 
momentum bias required in the above scheme and the effects of more 
than one momentum wheel, etc., are reserved for later sections. 

Refer to the section titled "System Performance" and see Appendix C. 
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STABILITY CONSIDERATIONS 

The desired action of the control law just developed is, 
of course, to reduce the angular velocity of the satellite to 
zero. Such a condition is not, properly, an equilibrium point 
for our system, since no provision has been made, as of yet, for 
a null state for the magnets. They are always on, having one 
polarity or the other, and the satellite always has some angular 
acceleration, hence the system possesses no equilibrium points. 

It would be easy however, to insert deadzones in the magn*- u 
controllers. Provided these deadzones are large enough r_pi. ^nt 
magnet activation due to the relatively small term dSPA (B) i,. 

dt 

(7) (which the satellite experiences due to its orbital_motion and 
the rotation of the Earth), then, the condition u> = 0, B = any, is 
an equilibrium point. If this sort of provision is not made, then 
it would seem that the satellite should behave in some sort of 
quasi-oscillatory manner, at small angular rates. That this con- 
dition would, once attained, persist, seems quite plausible from 
the analysis above. It is not, however, a simple matter to talk 
about equilibrium points and stability. The equations describing 
this system are quite intractable to ordinary stability analysis. 
For one thing, the control law is binary. Even if we replace it 
with a linear law (the first term of its describing function) the 
resulting system equations consist of six first order nonlinear 
differential equatioi ; three involve quadratic combinations, the 
other three cubic combinations of the state variables (w^, 102 » w 3, 
Bj_» B2 t B3) . 

We should note here that if the system has either a binary 

control law with deadzones or a linear control law (that is, M = 

• 

-KB), then there are an infinity of proper equilibrium points. 

These points belong to one of four distinct classes. However, in 
order to spare the reader, a discussion of these points and their 
stability is relegated to Appendix D. Suffice it to say that 
some of these points (other than w = 0) can be shown, by analysis, 
to be neutrally stable. Further, simulation results show that in 
some domain some of the points are stable. 

There are other problems too. We are particularly interested 
in investigating the existence of any so called "psuedo equilibrium 
points" (PEP's for brevity). Since dSPACE (§) ^ 0, and there are 

various disturbance torques acting on our satellite, and the 
system equations are nonlinear, there may be points which are not 
properly equilibrium (stationary) points but near which the system 
might "hang up". Here we are talking principally about limit cycle 
behavior, but must also be concerned with various (as yet unknown) 
forcing functions. 

In view of these difficulties, it was decided that digital 
simulation was the most suitable tool fc * evaluating this system. 

A discussion of .the results of this simulation work is reserved 
for the next section and Appendix E. 
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Although the actual stability analysis is put off until 
appendices C and E, the conclusions and implications of those 
analyses are reviewed here. First, we note that, for a vehicle 
configured as described in the section "Simulated Satellite" 

(that is to say, the proposed vehicle) , the terminal condition 
(or behavior) corresponds to a stable equilibrium point. This 
is, of course, most desirable. We again note that the "upside 
down" orientation (that is, with the pitch axis flipped over) 
corresponds to an unstable equilibrium state; this too is 
soothing. Finally we note that, for our vehicle, a spin rate 
of up to 27.8 x 10 -3 radians/sec. about the pitch axis also may 
be a stable equilibrium state. 

This stability conclusion neglects the fact that the direc- 
tion of B in space changes as the satellite orbits the Earth. 

This turns out to be a saving grace. Nonetheless, for strong 
enough magnets, the satellite may "track" the field. During such 
a period it is possible that the despin rate may be reduced by as 
much as a factor of fifty. For this reason it might be advisable 
to insert a mode switch for disabling the despin system for an 
eighth orbit or so if_it seems as though this has occur ed. At 
the end of that time w and B would no longer be parallel and hence 
the despin system will again perform properly. 
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SYSTEM PERFORMANCE 

Having convinced the reader (hopefully) that our system works, 
we now attempt to answer the question - how well? As before, we 
shall first be concerned with the despin phase. 

As the system kinetic energy ^ w‘I*w tends towards zero (this 

was assured in the section titled "Control Theory"), so must u and 
the system momentum f. 5 (this, of course, excludes the wheels' 
momenta). in fact, the despin process may be looked on as a re- 
duction of the momentum term f.j to zero. In the very best case 
this requires a time At 0 pt given by 

Atopt - LIihI - |Pq| ( 10 ) 

T max MB 

Here P D = I*wi n it t ^ le initial momentum of the body other chan 
that due to the wheels. It is useful to consider an "efficiency 
factor" a, 

= A tppt (11) 

At actual 

where a of course lies between zero and one. a is a function of 
various parameters, such as wi n it, the inertia configuration of 
the satellite, and the altitude and inclination of_orbit. That a. 
is less than one is a reflection of the fact that M is not always 
perpendicular to B and that t does not always lie opposite 

For a given system onfiguration (I & M) and a given orbit 
and initial conditions (determining u, B) we should be able to de- 
termine a and hence the time for despin. This was the goal of much 
of our early simulation, in which various orbits, initial conditions 
and system configurations were evaluated. This simulation shows 
that a generally lies in. the range 0.55-0.80. For the worst case, 
a near equatorial orbit with initial spin along the field vector, 
a may be as poor as 0.15-0.20. 

Another way to look at the performance of this system is from 
an energy point of view. Consider, for simplicity, rotation about 
a single axis. 


T = 

|lu 2 


( 12 ) 

T = 

I(|)U = WT 


(13) 

• 

• « 

• 

Iu> = x = aMB 

( 0 <a<l) 

(14) 


The last of these equations says that to, and the angular momentum 
Iw, are always decreasing. The constancy of this rate of decrease 
depends on the constancy of a. We ara of course assuming that t 
has the proper sense. 
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The attitude acquisition part of tie problem is more difficult 
to analyze. If, for example h approaches 0 we will not acquire. 

But also if h approaches infinity we will acquire, but it will take 
infinite time. Also, if M approaches 0 or infinity we will not 
attain the proper attitude. In view of these considerations it would 
seem that there might be some ratio of momentum bias to magnet 
torque which would provide optimum performance. By this we mean 
both minimum time to acquire and minimum roll-yaw error after ac- 
quisition. Consider the parameter 

Tsys = h (15) 

MB 

which has the dimensions of time. It is rather difficult to say 
anything about this problem analytically, but the following is 
observed from simulation runs. 

We first note that there is a certain minimum bias required for 
terminal attitude acquisition. It is shown in Appendix C, eq. (CIO), 
that 

h 0 > n max(I 1 ,I 2 ) (CIO) 

for us, ho = h + 138 and fl = 2 wo (for the terminal condition). 

Thus (CIO) becomes 

h > 2u) 0 {max(Il,I 2 ) - I 3 } (16) 

For good performance this amount should be exceeded by a factor of 
four or more. Little is gained by increasing h beyond this point. 

If we now choose a particular orbit, the time and quality of ac- 
quisition (after despin has been completed) are functions of the 
magnet strength M. Acquisition time is also, of course, very much a 
function of attituae at the conclusion of despin. A (very) qualita- 
tive picture of this situation is presented in Figure 1. As pointed 
out previously, for very weak or very strong magnets we do not attain 
the desired attitude. The author cannot caution the reader too 
strongly against taking this curve too literally. Much more work 
would be needed to get a more exact picture of this phenomena. There 
are too many variables and too few points on the curve of Figure 1 
to let the author rest comfortably. As for the quality of roll-yaw 
acquisition, it is interesting to note that a curve qualitatively 
similar to Figure 1 may be used to describe this as well. Figure 2 
depicts this situation. Note that the minimum point is shifted to 
about 2000 secs. (This appears to be significant) and the curve is 
somewhat flatter. The best that we can do appears to be about 2°. 
Again, the author cannot over-emphesize the tentative nature of these 
findings. Considerable work remains to be done in this area. 

For a well proportioned system, acquisition times are of the 
order of one orbit. This being the case, it is expected that 
optimum performance will generally be defined in terms of quality 
of roll-yaw acquisition rather than in terms of minimum time to 
acquire. 
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RECOMMENDATIONS FOR FURTHER WORK 

The author feels that the woj.1. done to date definitely demon- 
strates the feasibility and desirability of this acquisition con- 
trol system. Nevertheless, work remains to be done. 

There is the matter of system stability. While it seems 
apparent that there are no stability problems which might interfere 
with the proper operation of the control system, nonetheless it 
would be good if we obtained some further analytical results. In 
particular, an analysis including the effects of the change in B 
due to orbital motion would be interesting. 

As noted, we have not included the effects of disturbance 
torques on the performance of this system. Their effects on both 
limit cycle behavior and terminal attitude orientation should be 
investigated. Their neglect has been justified up to this point 
by the fact that the magnet torques are at least an order of mag- 
nitude greater than the disturbances. 

Another area which needs refinement is the magnetic field model 
used in our simulation. So far, only the simple tilted dipole has 
been used. The effect of a more accurate field model on terminal 
attitude orientation needs to be investigated. 

More investigation into the effects of deadbands, hysteresis 
loops, and inhibit controls placed in the magnet control loops would 
be desirable. There is also the matter of crosstalk between magnets 
and magnetometers. We must be sure that no unstable loop is set up 
in which switching magnets trigger magnetometers which in turn cause 
the magnets to switch again. 

One other point which will bear further investigation is the 
potential of using this acquisition system as a complete attitude 
pointing system. Used with a pitch reaction wheel scanner and simple 
pitch control loop we would have complete three axis control. Some 
studies of pointing accuracies versus various system parameters would 
be useful in investigating this possibility. 
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This appendix consists of a listing of the digital computer 
program used to simulate control system/satellite behavior. The 
program is quite similar to the one discussed in reference 1, and 
the reader is referred there for a more detailed discussion, flow- 
charts, etc. Its operation may be briefly summarized as follows. 

The main program performs the bookkeeping functions of input/ 
output and maintains overall control. What we are doing here is, 
essentially, integrating equation (1) , specialized for the vehicle 
described under "Simulated Satellite". This is done through sub- 
routines WDOT and RK4 . An additional problem we must deal with 
here is thhu the applied torques are functions of the vehicle's 
attitude. Namely, they are functions of the Earth's magnetic field 
components as measured in the vehicle. For this reason we must 
keep track of the vehicle's attitude. This is implemented through 
the use of Euler Homogeneous parameters in subroutines EHPQT, NORM, 
MAT1, and MULVEC. The Earth's magnetic field is first obtained 
in trajectory axis components (subroutine BFLDS1) , and then trans- 
ferred into body axis components, making use of the direction 
cosine matrix (subroutine EHPA) . Subroutine TORQUE then computes 
tne torques arising from the interaction between this field and 
the on board magnets. Please note that the subroutine displayed 
contains neither deadzones nor hysteresis loops, although some 
subroutines incorporating them have been used at times. 


t 


B-l 


T 

-J 


3 

3 


* 


‘1 

-j 

“1 

i 

2 

'i 

3 

3 

3 

3 

3 

3 

3 

3 

3 

i 

3 


r -4 CM ID 

o o o 
o o o 
o o c> 

CM CM CM 
> > > 
< < < 
o o o 


1 “ 

r *—# 

■> 

of _l ft 
UJ O ft 
_l > ft 
3 Z ft 
UI •— • > 

I * 
to y) -'i 
iu S * 
toui ft 
3;-J 

:CQ 

iu,o 

O.o: 

< Q. 

*j 

O'ai 
<> 
CL_I 


|» 

Of to 


;r 


< 1 *- 
ot! 
o oc 

UJ LLI 

*-,o 

z,of 
WiO o 

i 

o'z H 

< 

I-, H 

< l/> Q 
X af of 
UJ UJ 
Z H O 
»->;uj uj 
sf ,x ►- 

l< IL 
uj or. •-> 
to < o 
O 0 - lu 
cl! a 
of to to 

o.;o 

UJ 

z 

UJ 

o 

3 
X 

o 


< 

a: 

UJ 

z 

Ul 

o 

•c- 


LU 


° i 

y o u 


i 


in o 
o o 
o o 
o o 

CM (\l 
> > 
< < 
o a 

* ,* 

* ft 

* > 
*,# 

*;* 
ft * 
* 

* 

*•! • 

* o 

* H 
<- UJ 

* , 

ft - 

*- y> 

* oc 

•ir UJ 

«• :h 

* UJ 

* X 

*;< 

* ot 

* ,< 

ft cl 

*1 
*! - 
ft ,to 
* 

* ,< 
ft.H 
ft UJ 

ft ,o 

*;z 
*• •-* 
*i 

v 'j 

ft u. 
* 

*j 

* 1 • 

*,2 

* O 

* ;h 

* UJ 

* -j 

* UJ 

ft ;*r 
ft to 

* 

*- 
ft 
ft 
ft 
* 

» 

•li- 
ft 
ft 
* 

■ft 
ft 
ft 
ft 
ft 
ft 


H CO Ov’O ~4 

o o o y 
o o o o o 
o o o,o o 

CM (M CMCM (M 
> > > > > 
< < <!< < 
O O 0:0 O 
ft 

I* 

* 

* 

* 
ft 
ift 
* 

Ift 
ift 
:* 

► # 

< * 

Q-ift 
Xift 
UJift 
I ft 
► ft 

X ft 
Cfift 

ft 

*,* 
r-l ft 


ty> 

a c 

► UJ 

XX 
cc t- 
op 

to z 

UJ O 

zx 
— .< 
H; 

3 ► 
O, 0 l 
of fit 
CC UJ 
3 H 

wz 

to; 
uj o 
o z 
3 i< 

d ; ► 


o 

UJ;Z 

o < 

d'UJ 

_J 

e; - 
x 


o oo o 


CT 

a. 

x 

U! 

ft 

ft 

■li- 

ft 


o o 


(M ro -J- in -o H co 

H H H H,rJ H rJ 

o o o o o o o 
o o o o o o o 

CM CM CM CM CM CM CM 

> > > >,> > > 
< < < << < < 
O C O 0.0 o o 


CO 

to 

>c 

< 

r 

m 

•» 

CO 

< 


o 

o 

r- 

in 

o 


ft-l 

o 


ro.ru 
o o 

*'iu 

rvJiy 

0 qo 
l> w 

-I'i 

01 o 
O 


^ P — UJ 
13 03 


Jl 


m o x ,o 

>J- ~ ► X c 

► a h o o 
>J- 03,0 CJ "V 

I — ►>- 

a? m!m — uj 

— — m'x 

► i — ,3 — ■ 

mm ► i — ► 
•-> : v. 

z y -vlo 

► *-«—* l/) e 

fNI j a H O 
O V. 3 V 
Z * O uj 

— 0 £ 3 X 

H 'N, >— ■ 

•- sj- Z Z , H 

■ii o o, 

_l _l X X _J 
< < X X < 
jj uj O O ui 
o; at o o or 


ru m '*■ in >o h cca>o 


O' o •— * ' fNI 
.-i rv cnj csj 
o o o'o 
o o o o 
nj nj rjp 
> > > > 

< < <;< 
O O 0.0 


m y 
nj r\j 
o o 
o o 

CM CM 
> > 
< < 
o o 


ar 
o 

H 
O 

if i 

► cl! 
to UJ , 

< I— , 


CD 


m 


CM 


to 

CO’ 

•» ' 

Hi 

Z| 

M 

cc! 
a. ' 


o o 

O' O' 


in in in 


o o 
< < 
UJ UJ 
Of of 


m 

ft 

« 


< 

—> * 
;< f-» 
1 *• 0 
Of O 

X 

H H 
Cf Of 

.< o 
UJ to 
Z ► ft 

3 ,z m 
ac —i ~o 

H J CM 
••OH 
H Z m 
G *-• r-C 

m 

0-0 0 
O' 'O' o 
► ► n 
m uo o 

o 

a o — 

< < Of 

Uj U! UJ 

or cl' a 


in <o H 

CM C\J CM 
O O O 

o o o 

CM CM CM 
> > > 
< < < 
o o o 

I 

I 


Report 
Page B 


No, 90,345 

f 2 


co a o 
cm cm m 
o o o 
,o o o 

CM CM CM 
> > > 
:< < < 
O O O 


O 

O 

M 

\<X 

& 


\C£ 

.X 

,H 

Of 

< 


Of 
O 
H 
O 
< 

U- 
•» 

< 

•— Q.UJ 
CO UJ, ► 
► H Z 
m m ji— < 

*-< <D _J 

o 

CM HZ 

— z ►- 

► 

of :h 
-« 0 . 0 


cm m 


o CM 
H f — 4 H 
H H H 
•> •• ► 

vO •O.^ 

LU UJ UJ 


a or of 

3 3 3 


m r-l CM 


i l 
m f-t 


ii ti 
cm m 
i_> o 


cm miy in 
m m m’m m 
o o oo o 
O O 0,0 o 
CM CM cm CM CM 
> > > > > 
< < < < < 
O o op o 


UJ 

o| 

M 


H! 


i 


■c h cc O' O r-i cm m 
m m m m|^- -f- m- 
o o o oo o o o 
o o o oo o o o 
cmcmcmcmcmcmcmcm 
>>>>>>>> 
< < < < < < < < 
O O O 0.0 op o, 
i 1 


UJ 

of 

UJ 

3 

to 

Cl 

X 


uj a. — 


o 0 s 


o 

CL' O' 

H cm m 3 ► 
i— i w « in 


u. 

* 

1-4 

o 

o 


(O 

io 

I *— < 

•OC 

is 

O 

o’ 

o 

m 


*— (D 
H 1 ► 
CD < 
•■')- 
H H U) 
O -J X 
ft <,H 

o y ► 
o Z Of 




i 

— UI 


r-l r-l H Z C 

II II ' II < < 

-fr in «o iu 

o o o o of 

<! 

UI 
of 


UI < 

HiX 

— ifif 

cf io 


m|m M-'in o h co|o- oh 
- 1 - 1 (M cm 


H H J*H H H H H H 


i 

.1 I 


cm m 
CM CM 


CD 
- II 
O — 

II r-l 

H ^ 

to >: 

z x 
o o 
o o 


Q. 
UJ 
' H 
* tO 
II CD 
— , II 
CM CL 
— UI 

X H 

X to 

o ae 

u — 


Of J.H O 

uj o of O 
Q. Z < *-* ! 
V.*-*UJOf‘ 
H |0 w w U) 
z cHNaj 
H O to'to >N 1 
JCOCQN 


m —i 
ii u 
< a 


_j m 
u. co 
CO CM 


1 - 4 ! UJ 

Cf i 

CL 


O' 


y n h cc o o r-i cm m 
CM CM CM CM CM CM 'CD m m m, 


-I -J o 
X -I -I II 
H < < O 
O O O 3 , 


7 >, '• ii. oil' mid 


I Ilf Cw mscu 


I' _ 





in \o f- oo 
c o o o o 

o o o o o 

CM <M CM CM <M 

> > > > > 
< < < <r < 
0 0 0 0 0 


O' O rt (VI 
N * in in'iA 

O O o o 
0 0 0,0 
CM CM CM CM 
> > > > 
< <r <>< 
O Q OO 


O 

< 

S' 
a 
in - 
** a 
o- 
* £ 
t- a: 
oo 
n z 
u. U 

< _i + 

- _j L- 
*- < ii 
o o' 


in 

r- 


o 

o 


UJ 

z 

o 


m 

« 

m 

3 

* 

m 

M 

♦ 

CM 

2 

* 

CM 

3 

* 

CM 


I 

cn 

in in 
o o 
o o 

<M (M 
> > 
< < 
o o 


m 

CM 

< 

► 

m 


- > 
mo o 

W CC 

UJ UJ 


in >o r- co O' 
in m in m in 
o o o o o 
o o o o o 

CM CM CM cm CM 
>>>,>> 
< < < < < 
a o oo o 

i 


'<* 

I 


Re 

Pci 


porjb 
gf 


Nd 

fc-3 


9 '0 3 4 5 


z ~ 


a,o 

•"■i ■ 

* A 
y-' i 
z o 

M '«£ 

a; X 

o-:a 

•-< ii 

i x-ac 

->C 

'-'at 

U. DC 


3 Q 

* ,< 

r-l’UJ 

— X 
II - 
>•! ii 
o^o 

lu'lu 

Z ; X U. 


M 

CC 

2 


O 

< 

Ul 

X 


o 

UJ 

9 

o 

in 

* 

o 

in 

\ 

o 

< 

UJ 

X 


M' m 

cn fc ( 


>0 n> 


oo cn 


O rH 


m cn 'cn ci'« 4 - ^ 


UJ 

•» 

aC 

O 

DC 

DC 

UJ 

•» 

CD 

i® 


LU 

x 


s 


* 
Ul 
2 
+ 
cc 
X ■ 
I— l 
DC 


CD 


CM UJ 

oo X 

o y~ 

_» ~r 
U- CD 
CD 00 

o 

-J _J 
_l U- 
< CD 
O 


< 

e O 

;cr 

UJ — 

0 fM 


■< 


< ; a 
h-;lu 
lu,h 
X:m 
J- CD 
C— * 


o o a 
o + — 
ui < 
X Cl 
O —.X 
I- UJ 


UJ II 


UJ _l 

— x;_i 
X —;< 
2 I- O 


m 

n* 


p 

<A 

H LU 
Ul H- 
X !oo 
t— ,co 
ii — ■ 

y- 

ur— 

x u 


s 

o' 

*i\ 

CO ‘ 
C Q ! 


[CM 

>4- < 


cn in!o r- 
^ •4" 


UJ 

z> 

o 

DC 

o 


<; s 

u o 


O —i 

CM 


-J- in 

<0 

f- 

CO O' 

o 

1 

r-l CM 

cn 

«4 

in 

1 i 

mO H CO O' O 

vO 

MO 

O 

o:o 

*0 

o 

o 

mo mo n- 


K 

n- 

r~ 

r— r- n- 

r-« ® 

o o 

o o o 

o 

o o 

o o o o o 

o 

o 

o 

o o o 

o o 

o 

o 

o 

o 

o 

o 

o o 

o o o 

o o 

o 

o 

o 

o o o 

o o 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM CM CM 

CM CM 

CM 

CM 

CM 

t\ CM CM 

CM IM 

> 

> 

> 

> > 

> > 

> 

> > 

> 

> > 

> 

> 

> > > > 

> > 

< < 

< < 

< 

< 

< < 

< < < 

< < < 

< 

< 

< < <r 

< < 

o o o o 

o o o 

o 

Q O 

o 

o o 

Q 

a 

o 

O Q Q Q O 






GO 





1 


o 

LO 

<M Z 










m 


- 


z 

cc 

<« UJ 

in 









ii 


CM 


o 

UJ 

O Ui u 

in 









- CC 


0 


joe 

H- 

rH X - 

- LU 








•> 

CM O 


o 



UJ 

U H- 

fl -I 









L H* 


\r-4 


j 1 

X 

•. iu 

1 - z 






H 


II 

O O 


*',u- 


! i 

< 

- LL D 

— z> 






CD' 


rH < 


o! ► 


*■* 

cc 

II O 3 CD 







— - 

H 

U U 


* - 


loo 

< 

>- 

cc x 







« 


•» 


o> 


c 

d'O UJ 

O LU 





m 

I2C < 

- X 


HOC 


z 


Z Q 1— 

y~ 





0 


< 

* — i 

II h- 


U LU 


o 

t0'< ID —I UO 00 





r- 


>- 

i- 

— O 


-■> 


l_> 


l- < 

*— > 





U- 

CO 

o: 

cn z 


• UJ 


Ul 

O Z — ' 

X VO 





(NJ 

< 

v k 

UJ 

— Ul 


>1 


o 

UJ 

O 10 h 1- 





•o 

0.0 z 

< cc 


CC O 


j 

z 

— z 

oO 





ur\ 

at 

UJ 


h- 


UJ 1 — 1 



LU 

>— O CO 

cc 1 





•i 

z 

Z 


h- oo 


>IU 


1 

o 

<5 — * OO 

OX 





in o 


»■ 

DC 


LU 


m 

c 

z o u 1 





rW 



XUJh 


UL 


mt 

X 

^ CO - 

X 





UJ 


_1 

O 

Z UJ 


\ 


:cj 

o 

-I X •* 

o, 





cvi 


_l 

»0 

z 


z o 


I ® 

X 

O Z X Q z 







c 


- o 


f—t f—f 


'in 


Z — Q 

M M 





*> 

a: s 

► < 


cc y 


i 11 - 

DC 

h £ in 

a: 








V. 

X X 


Cl UJ 


! « 

L'f 

» ► LU UJ 





X 

CM 


•» 

o * 


- z 


'» 

. , 

► I \ Q. 0^ 





r-l < 


«■ 

•o - 


- o 


i •! 

o 

X J- X 

• < 






O 


UJ 

- X 


x,< 


LU 

Ul 

m oc - 

•k • 






UJ 


QC 

X o 


O X 


A 


CM < - 

X oo 





n 

X v 

< 

v. -n 


MO i 


*-4 

X 

*> UJ IA k 





o 

o o 


•. » 


-Ul 


;«n 

1- 

X LU 

CM *-> 





u. 


of in 

CM V. 


'v.H- 



Hi 

X U. U' 

- z 





>— 


iu ac 

r. N 


>*•< 


'Q. 

2 

► O CC 

X Z5 





m 


Z UJ o •• 


- -I 


|uj 


- o \| 





•* 


IU 1- 

r-< no 


- P 


V 

Z 

oo Z UI 

-— 1 








Ul 

UL 


ui;o 


00 

a 

z o a 

- — J 





X 



2 : 

- O 


ccU 


1 

r-r 

a — 

r < 





*-4 



< 

•* r-l 


<i< 


;Z 

I— 

►H 

X - 





#• 

< 


a: 

II U 




,c 

<r 

1- < 

X| - 





•n o 


< 

rs r* 


00- 


M 

cc 

< K* 

|X 






UJ; 

CL 

CM - 


-J - 


A 

o 

a o - 






o 




^ II 


<;x 


< 

IU 

►- CC CM 

CM 





-i o 


X 

< LO 


> in 

m 

CC 

*— 

U - 0 

- - 


CL 



UJ 



UJ 

—1 < 


CC CM 

k- 

o 

z 

r- < — O 

- A 


o 



w 



h- 

1 — ' 


Ul 1 — 


IU 

- 

O X r-l 

O v. 


H- 



CO 


CC 

uo 

cc co 


1— v 

CC 

1- 


ui m u. 

; M 


oo 



► 


o >- 

Ul 


z:n 

CL Z 

cc 

CL CM - 

O - 





X 

OJ OO 

Z X 


- 

o 

*— < 

UJ 

LO - - 

r-l - 


— 



<\] 

•QC 

«• 

-H Z) 


- * 


• 

a 

S II 

u. j C 


z 




UJ 

UJ 

•» 

- h- 


- m. 

u 

•» 

cc 

k v.;— 

- OO 


z> 




z 


X 

- z 


x A 

o 

X 

r-> 

— i - J 

- o 


DC 



o 



o 

x u: 


O UJ 


o 


a - < 

II UJ 


1 — 



e 



in 

O X 


uo 

OO 


X 

(C IS) r— 

— m 




CO 

nO 




MO O 


-oo 

UJ 

r« 

v— 

o un- 

u»j 


0 


r* 

LL 

• 



- X 


x! 

UJ 

S 

cc 

• LU M 

w . 


H- 


O 

•> 

► 



s . 


Ul 

cc 


p 

- CC z 

z , 


O 




H 


•m 

s •« 


X X O 


o 

X o •- 

< • 


V 


UL 

X 

X 


s 

‘ X 


\ r-» 

UJiN*. 

u. 

OJJ- 

H- — 

M- 

o- 

LA 

w—i 

r-im 

N 

CM O 

(M 


o 

,*s. 


in D in 

m o 

X 

'jj 

in 

w« 

W 

W 35 

•-* 

f> MO 

c 

•— * 1 



< 

- UJ 

—• 1 o 

OC X 



1- 

H 


h~ 

o ► 

O 

Hi 


:i— 

1— 

■V Q Q o 


*— i 

O «r 

< < 


< 

r-l "T 

-i < 


;< 

K 

- a 


-J H- 

H 

* 

* 

X 


X 

u. -V 

u- 

X - 

m 


Z) 

x - Z 

-J u. 


u o 

M O 


CO cn 

d >4 - 


DC 

o 

u 

o 

O' 


o rj’cM m 
m in in in 


a: 

o 

u 

o 

o 


ac 

o 

U CM 

!oa 


>- in 
m in 


° 
u. i 

O rU CM 

O. 


«Ol 
Ul 1 


CO ,C 0 

o o 

o o 

CM CM 

>;> 
< < 
0:0 


CC 

O 

U 

CO c 


UJ 

00 


r- 

in 


cc x 
jo i 
;u id 

N O H 

12 


DO 

in 


v • S> 

- cc 

- v- 1 

cm cn <j- in of' 


o 

z 

Ul 


O' 

in 


»T cv t 



r <0 «J- 
co co 
o o 
o o 

CM CM 
> > 
< < 
o o 


V. o 
co oc 
o o 
o o 

CM CM 
> > 
< < 
O D 


jC 

IU 

z 

iO > 

Z II 

'►-4 

o o 
•Z -I 

0 o 

CL — 

;co> 

;UJ » 
!C£ < 
loo 
lO 
JO « 
- 

<■ 

1 

IX O 
•— l/l 

c c 

H CL 

|< 

£ nO 

* 
Z UJ 
O H 
I »— • UJ 


I 


< < 

h CC 

o < 

00 a 

ui no 
xo 

H O 
UJ 

;no z 

— Z UJ 

< «'a: O 

►MDO 

Q_ »*, | — Tl 
-rtUJO 

< — CO X 

a <. 

X UJ co 
uj •> z UJ 


r- cc O' o 

co co 00 O' 

o o o o 
o o o o 

CM CM CM CM 
> > > > 
< < < < 
o o o o 


CMlcO 

O' O'lO' O' 
o o o o 
O 0,0 o 

CM CM. CM CM 
> > > > 
< < < < 
c oo o 


co 

!<m cl 

CO + 

A *- 
j I CM 
:CM - 
cm a 

Q. ji 

I 1 ~ 

«M r-M 


•J" «fr 

a o. 

■K- * | 

CM co! 


IT O H CO ' O' O 
O' O' O' O'lO' O O 

O o O O O rrt r-C 

o o O oo o o 

CM CM CM CM 'CM CM CM 
> > > >!> > > 
<<<<,<<< 
O a O 0;0 O O 

* 


f^epd 

ag4 


rt 

B-l 


(No. I 90 

■A 


345 


'*•!<• mt 


a,CL 
*• j* 

3? 


+:+ 

cn‘n 

CLCL 

* j« 


ui^hDaaa 

Z"3UJ'( v * 

— Q- o ! 

i— cc uj ; »— * .m n 

O M- I 

O -if wi-aaa 
ac _t — li ii ii 
CO <t X O N CM CM 
3 UJ H I — J CM cn 


r-* .M CO jo. CL ]0 l 

<x!cm * U 
* o s 
— a CM CM 
>* II II II 


CL CM CL A 

I I'-' 

^ Cl 
cm' I 

— CM 

a co 

iff 
3 H 

a a 
-‘ + 

# ;cm 

o 

CM CL 

II ; II 


2 L. 

tt. 

# j* 

£>! o 
CM <M 
ii ; ii 


CL 

I 

r 0 

CL 

»M 


i/i a. 


* 

:* 

* 


o o 


O 

O O' 


cl a CL 


CL 

n 1 •> ► 

NHH 
■J — — 
CL <T < 


CM]Cf| 


cm m 


CM 


CM ,<M CM CO cn 

< < < <'< < 


CM cAm - in O H;CO O' o f-l'CM CO 

o o o o 'O o No -oh h h h 


CM 

r*C 

CL 

I 

CM 

CM 

CL 

I 

CM 

CO 

a 

+ 

CM 

>J- 

CL 

II 


CM CO M" IO|o C* CO O' O r-c (M CO M- to o H OO O' 

O O O O O O O O r-c r-t‘r-c rt H rrt 'r-t r-C^rM rrt 

rH r-l r-4 r*^ H H H H »-t ‘f-4 H H H H r-4 1 « — < r-4 

OOOOOOOOO OO O O O o o'o o 

■ CM «M CM CM CM CM CM CM CM CM CM CM <M CM CM Cm'cM CM 1 

> > > >,> > > > > > > > > > > >;> >; 

< << < < < < < < < < < < < < < < < 

aooojoaoooaoooooo.oa! 


H* 

UJ 

L0 


UJ 

> 

>— I 

o 

>- 
z 
< 

CC 

o 
u. >rt 

UJ c&| 

LO O 

O no 
CL 

I/O CL 

z 

< UJ 
cc X 
h h 

O I 
I 

X I 
w LO 
OC CL 
I- UI 


I 


UI < 
X CC 
< 
a 

CL 

3 no 
ZJ 

LO o 
UJ 

M- UJ z 
► «/> UJ 

-t o 

UJ O 
HZ Z 
Or O 

itr x 

— o no 
<• oi * 

•— 00 CC 

cl'o uj 


•J" — CM co 

*»— ■ 1 w» W 

a i a a 


ii ii 


— CL 
CL I 
II II 


Z' 
co cC' 

- Z3 

m h o 
— uj z 
< cc uj 


r> 

o 

cc 

CO 

,r> 

no 


-j — ^ 

Nfj ro •-« CM CO rrt CM f*0 
•«-l/)Uj!r- •• . ► . ► 

_J»— rtrtrtNNN 


jCM- 
!— r-l 

CL — 

I CL 

II II 


— .-j 'cm CO 


X U. 


j- ICO C 
H H H 


H CO 
H H 


O O 


H J- .1- H 

a o o o 


H H 

a o 


N — 

— a 
a i 
ii ii 


A Ol 

l l 


rt n!cO rt ,cm co 


O' O .rt CM 


co c<0 JrO ^ jM - 

H H 

o o 


a o 


>0 M- [in M 5 


O O 


r- co 


O' o 


H oo co co co aojoo co ao oo oj O' 


OC 

o 

a 
uj z 
oc uj 


CM, 

O' O ' 1 


I 


) 



I 


Report Nc . 9p345 
Page B-6 


cm m sf in'o r- cc o' o — • cnj 
m m m.inp so >0 

*-4 »— 4 H H *— 4 H H »-4 H 

o o o o o o o o o o o 

> > > > > > > > > >;> 

< < < <.< < < < < <,< 


ooo 

ooo 


* 


* 

A 

«■ 

>■ 

* 0 

1- o 

* H- 

ri 1— 

* or 

o o 


o + ! 

* > 

— 1 h~ 1 

•H- < 

IU 

0£ 

> o! 

* DC 

h- 

* j< 

cc \ 

* 

< H 

* o 

_J 

* z 

O uj: 

* < 

O X 

* 

z — 

* •- 

<1— 0 

* ]•— 1 

UJ 

<-» -r «e- 

* *- 


m ,-4- in v0 
'O |<© 'C 'O 

H H H •— I 

oo o o 

IM (M CM nj 
> > > > 
< < < < 
O O Q O 


f- cC O' O *-< CM in vO f- CO 

>0 >C >0 N S f* N f— f*- 


OlO o o o 

CM CM CM CM CM 
>,>>>> 
<!< < < < 
O p Q O Q 


r- ,r- — r'^ — r 

r4 r*H jr-4 r-4 r-4 H 

o o o o o o o 

CM CM <\i CM CM (M (M 
>>>>>>> 
<<;<<<<< 
O 0:0 O Q O Q 


O' o «-i tv; fr> in Iso eo O' o 

r-cococooooocojcocococoO' 

r-4 *—4 H H H *-4 H H H 

opoooooooooo 

>>>>>>>>>>>> 

<.<<<<<<<<<<< 

OOOOOQOOOOOQ 

i ! ! 1 


cm m 
, O' O' cr 
^ .-c *-c . 
ooo 

CM CM 00 ’ 

;> > > 
< < < ■ 
O O O c 


O iu O uj 
0. 0 OC;X 
► Z u. o 

— UJ It/) 



X < n't 

~ IU 

-OMZ 

< . *■ o •-* 

— ► t- 

a. >4- — O 

— m o 

I Vo — cC 

yOrico 

0:103 

s m 
2 >t 

1 a # * t/> 

7" _l _l — • 
' X < < X 
O UJ ID h 
IO0tC£ » 

I * 
* 

s # 


o 

— 0. 
h •* *** 

O 3 H 
•* •* c 
cl h ► 
-C3fO 

K w | r. 

O r-C V-!r-t 

o. f- O ’ M 

x < o,~ 

W X X ! 

i© 


< < < c 
0000 


r-l Q 

O 0. 

* ft 

h 1 I— 

o | o 

ft | ft 

in in 

»! 0 

o p o 
+ , •> + 


— ;c — 

moo 


a. m o 

* o •» 

1 — — 0 . 1 — 

a *-, ► *•* o 

k- o m co 1 
in •> o 
0 o )- •-! 

*J- O - C 
»■ + I— —* } •> 4 . 

r-l ~ O r-l H r-l — 

11 Mat*C|ii — 

*-« ■»» X < Q — — 

a uj £ 3 i X 

*-•: 11 ] p 11 

<m ~ j j j ( n- 

H J J J' M 

o—<<<p~ 
o o o o 00 m 


o ~ 

cl ! <• 

«■ ! o 
H,-0' 
o t- 

a mp 

; •» •» O 

Q h' » 

^ — • o m 

► + i- — ■ — 

H ■ 1 ■ O 1 — l ,h“ 

11 w a h o 

l-H X < O 
Cl UJ X X 

•— * 11 : ! 

m — _j _j -J 

•— j 1 1 

O w < < < 

o o o o o 



:>t - \',h ef\ 

,0 * w »-» «— P 

»• h a a o. o. 

IfllflO iH # -if 

*-+0000 
,f- r-c — X x 3S X 
O II — + I l ; + 
a 1— 1 - 11 11 n , 11 
x x — — 

O # HPOCD,^ 

-J M - — — ,w 

_i -.0000 

<0-0000 
o c 'x a a a ,a 


lo o p o p 


O' O r-l CM 
•*< CM CM CM 


| | | I | I * I 

K) 4- m,«0 h CS O' O rC (Mio 4- Ift -0 N C6 O' O HW (O^ift.^heOCC^N 
cMCMCM(MCMCMfM'nr".f'. cnmnlcncnco My >#■< m- >4- -finmin 

H r-l H 'rl H .»! rl —I 1—1 ^ f-l H H |rl »— C rU rl rJ rl jt— I r— < r -4 rl 'rl r-l |H t— ^ e— - rJ H 

I i ! i 


CL Z 
O II CC 

in — o 

0 *-< I- c 

1 O — UJ z 

o D a c: i: 

X 

O 

in 

•if 

* 

CC 

o ! 

cc m < m l 

o; in m in i. 

U! r-C r-t r-l r 

■> i 

* I 

* ■ I 


? < ft ni aiis5>*4 


u.r » ON 


I 




I 

! 

fA vj- 

ff> m 
fM fM 

o o 

fM CM 
> > 
< < 

0 o 

1 


ia 'Q 
ro m fi 
CM fM fM 

o o o 

fM fM CM 

> > > 

< <,« 

O 0,0 


to 

O 

o 

k-l 

> 

UJ 

of 

& 

O! 

to| 

> 


't/I 

2 

■ ► 

to 


oc 

UJ LL' 

p X 
a o; 
a: >- 
;o *— 
p 3 : 

! to 

UJ 

2 H 

— ■ UJ 

H 2 
O O 
O < 
Of 2 
CO 

o _J 
CO < 
UJ 

o 


0 

CO 


to 

H- 

O UJ 
O VO 
UU 
£0 of 
a 
o 
I- o 
2 

2 -< 
O Of 
— < 
to a 
2 
2 o 
*-< o 

Ul > 
p CO 

'to Q 

O UJ 
|CL H 

[O 2 
! of 
!to uj 
uj j- 
•- UJ 
K O 


to 


a to 
A H- 
o 
h- o 

UJ 

2 <r. 
o 
< 

2 


to 


I 


UJ 

to to 


0 ;'_> o 


o 

•0* 


CO o 
fO ro 
fM fM 

o o 

fM <M 
> > 
< < 
o o 


O'rt fMifO 
•4- 4- sT ;4’ 
fM <N fM fM 

O O 0:0 

fM CM <M,(M 

> > >;> 
< < < < 
o o op 


UJ 

*-t 

u. 


UJ 

2 

o 

< 


u. 

o 

to 


o 

e 

O 

•M 

o 

_J 

o 

m 

CO 


O 

0 

o 

•v 

o 

•o 

1 ™ 

CO 


o 

o 

o 

a 


CD 




h ce ;2 

2 

UJ 
2 
C 
O. 

2 
O 
o 


to, •• 
to — 
Op 
O'" 
O CO 


to 


I 


II ffl 


o 
a 

cf 4 - »-f cm 

O 4c — ^ 

t- -j ca cc 
;< 11 n 

»• LLI «-< CM 

2, of eo ce 

uj! 

X 


4 'm o>* 

4 4" 4^ 4^ 
fM fM (M <M 
OOOO 
fM CM <M fM 
> > > > 
< < < < 
OOOO 


CC O' Q ,H fM fA 4 


1' o 


4. 4 * in a tn m tn; ■ in 

fM fM f\J (N M fM fM.fM fM 1 
OO OL JO 0|0 O 
fM fM <\| CM M fM AJ fM CM 
>>>>>>>>> 
<■< <•■*<< <;< <; 
000000 op o 


ji Report No, 
l| Page B -7 


tA sO f» CO 


o!o o 

UJ’UJ UJ 
2*2 2 
*H } CM f«1 

op o 

CO Ci CD 


33 

— P 

CD|2 

II 


U U 


CMp 
CD, CD 

I I 
0,0 
-J _l 
OO 
CM (A 
CD CD 
It II 

o o 

UJ UJ 
2 2 
fM CO 

o o 

CD CO 


p fM.fO 4|lA <0 V* 03 
qs O' O O')© 4 O' O' O' 


to to 
2 2 
o o 
o o 


CM CA, 

CD CC cD 
A * * 
n f*s cm 
2 2,2 
I 111 

fO r-tCM 
CO CO CD 


3 * # 

2 2] CM CA 
0 0 2 2 
—• •- It II 
to top ~ 
II II 
fM r, 


H CM HI 
CO CD CD 
II It II 
OOO 


2 2 


O' O 
O' o 
r->. CM 


rtNOQOO 
H(M(f| 
CD (D CC 


I 


*-< cm ro «t in <o 

000000 

fM fM fM CM <M CM 


2 1 

§ I 

!uu 2 

of UJ 


1 


r~‘co 

o o 

CM fM 


, 9 i 3 4 5 ; 


1 


O' O r-c fM fO 

O' O' O' O' (A o 000 

Hrlrf ff'w fM CM fM fM 
OOOOOOOOO 
CMfMfMfMfMCNifMCMCM 

> > > >;> > > > > 

< «i < <,< < <f <r,< 

OOOOOOOOO 


— to 
.CM — 
> CM 
*>, 


< f^l 

1 w w; 

o «-< 

UJ > 

> N 

=J - t > 

ip — ! « 


2 fn'r-e 

-l> 


'UJ ff|,* (A 
'2 ► 
•-« <, r-t 
f- to II 
,0 ^ UJ <-l 

0 # *-« 

Cf _J _J 

co < a 
D tu « O 
t/J ffi K- O 

1 P 

2 


O' 

tn 


f- co 

tn tA’ 


Is 

I* 

I-* 

r* 

1 ^^ 


fnU 

o •«!*- 

OHM 

O II ,> 

II Sf 11 2 1 

— p Of 

M H M 3 1 
w ;vi-p 
fM C CM UJ 2 
> O > Of UJ 

i 

i 


O CM fO p 

^ vQ 


i 



1 





B0345 


) 



•J 

I 


3 

i 

J 

i 



j 


i 

i 


) 





■*r 


i 


«j- un 
o o 

CM CM 

o o 

fM (SI 
> > 
< < 
Q O 

UJ 


I 


M3 f*|0 C O'OhNIO 

O 0,0 O — c — c 

CM CM CM fM CM CM 

O 0,0 o o o 

CM CM, CM CM CM <M 

> >,> > > > 

< <,< < < < 

O O O o o o 


H H'H H H 
CM CM' CM fM r\J 

o 0,0 o O 

eg (M Tsl fM <M 
> > > > > 
< < <t < < 
Q O P O O 


< 

QC 

< 

a. 

ILI 

i/i 

UJ 

IU 

1/1 

I 

I 


CD 

1C 

o 


V) 

z 

o 


i/1 

o 

<x 

ul 

-O' 
.to — «| 
* Ot ' 

,< < 

, •• > 

z 

cc 

-I o 

o u. 

z 

O 

!— ■ -j 

1-c Ul 
: (/) •-« 
O LL 


UL O 
<D M 


IU UJ 

z z 
— o 
>* < 
3 X 

o 

or ui 

to Z 
3 a 
u» 3 
»- 
UJ 

! tz 


o 


o 

rH 
-* 
f-4 
e 
(ft 
*s 
— » 
a 

•» 

z 


O 

z 


Z eft 
O ~ 
— * CO 


< 
* 
< 
* 
< 
# 
— « 
ft. 
* 
o 

e 

<4- 


O O 


UN 

>0 


£<s 

X UJ 
3 at 
o 
O 
o 


ftj 


il 

—4 

O.fM 

e o 

CM 
•I , e 
ST O 
OH 

z z 

CO Ul 


UJ IU 

p o 
p o 

!<t < 
ia at 
!* * 
;Z Z 


o o 
'z z 

O M M 

CO,— — 

O' Z L 1 

. c 
O U 1 o 
II II II 
Z,« -*iUJ 
O eft o'k 


Moa>.o 

r-i f-4 -4,CM 
CM CM CM , CM 

O O 0,0 
cm eg cm 'cm 
> > >!> 
< < <:< 
0.0 op 


1/1 
UJ 
UJ 
I— at 
y> O 
uu UJ 
uj o 
lot 

p O 

UJ V) 

Q ~J 

<c 


UJ 

a 

o 

z 

p 

z 


o 

z 


a 
z 

UJ 
o «— • 
VI O 

< z 

UJ 

sc o 

O ui 

at < 

UL 

UJ 
O X 

Ct h- 

< 

25 U. 
— Ot O 
CO o 
► ;u. uj 
< o 

H iK- 


O 


e* co O' o • tu 
so «o «o r- c- r- 

rl H H W H rl 

I 


3 

tlj • K 

X;co — 

I— IOC OiUJ 

•■'a zp 

at; Op 

x;z _J < 

at i j* 

< |UJ H* i— 

UJ 1_J * at 

— p — |x 
eg ;z X ,i- 
ui,< kiq : 

O fc,< 

J |uj < UJ — 
u. x uj : i o 
CD h — p — 

I X, * Oft 

> v> O' 
ct *- it ,o 
j-, ct — 

■ ' x ; II 
1-3 
ot 
< 

UJ 


z 

LJ 


S 


p u 


M >t 
/*■ 

f-J 


o 

o 

ii 

3 

O 


»*l (M (ft 
(M f\J (\J 
(M (\J CM 
O O O 
CNJ (M (M 
> > > 
< < < 
o o o 


<4- iftp 
CNJ (Nj'cv 

<M (\J ,r\j 

o oo 

CM CM <M 

> >,> 
< < < 
o op 


z 

i/1 

* 

3 

t/l 

» 

M 

VI 

+ 

Ol 

* I 


* ! 
z z 

O Ul 

z « 
to 3 
II O 
V) — | 
II (ft II 
3 — r-i 


>1 


Z 

0 
■* 
t-4 O 
i/1 LU 

1 o 
z o 
to < 

«■ at — 
3 •Jr > 
wi < — 
* J- V) 
« m O 
X o 
II 

II > 


o 

II 

CM 


Iflef oo>o 

I 


UN«C|C-CG0'O— «CM 

p* r- r* co ® co 


P- 00 
fM CM 
CM (M 
O O 
CM (M 
> > 
< < 
o o 


eport No. 
kge| B-d 


O' O — 1 1 CM 
fM (ft fft (ft 
fM fM CM, CM 
O O Op 
CM CM fM fM 
> > > > 
<< < < 
o o op 


> 
t/l I 
* 
cm;. 
o;> 
+ ;o 
> # 
O fM 

* o 
-c I 

o> 

— VI 

* ;* 


|UJ 

lx 

h- 


y> 


o o 
z,— 

CD ’«• 

* X 

o o 

o Z 
fNI CO 
II III 
ul — — 

II M fM t- h- 

> 

Ul CO CQ 


> 

Z 


I t 


,h- r- 

.t-4 w* 


Cft 

CO CC 


u. O 

(O r-i 

Z I 

Ot O 

3 .at 


ujiZ O 
ot iuj O 


IP gj NtO 
CC CC oo CO 

[—4 — * — I — 4 !r4 —4 


o 

z 

UJ 


f 


* ui ui ‘jjiumj 


1 1 I f i ON 




Report No. 90345 


APPENDIX C 


DUAL SPIN EQUILIBRIUM CONSIDERATIONS 

* 
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In this appendix we consider tiu ...x^ibrium and stability of 
a dual spin configuration. Although this work has been done before 
(e.g. reference 6), we derive the results here in a simpler way. 
These results were obtained independently by the author before his 
acquaintance with reference 6 and are included here for complete- 
ness. 


Consider a body with a single momentum bias wheel; the wheel 
spins at constant rate and carries angular momentum h = {0,0,h}. 
The body is spinning about the same axis as is the rotor; for 
definiteness this is taken as the " 3 " axis. Suppose we consider 
the torque free motion of the body in the "neighborhood" of an 


equilibrium point , namely 

(i>l = u>2 = 0 , u )3 = ft (Cl) 

In the neighborhood of this point 

• • 

u>2 = 6w^, = 6(1)3 (C 2 ,a) 

(i) 2 = b>2 = (C 2 ,b) 

&>3 = J 2 + 60J3, = 6(1)3 (C 2 ^,c) 

Now, the equations of. motion, adapted from eq. ( 1 ) are 
I3W3 w 2 ^ (I3 ~ I2 ) ) e 0 (C 3 ,a) 

I 2 w 2 - ^{h + fc> 3 (I3 - Ii)} = 0 (C 3 ,b) 

I3W3 + (1)3 u 2 f 1 .. - ^l) ~ 0 (C 3 ,c) 

Substituting eq. (C 2 ) into these, and dropping second order terms 
(involving products of the 6u's), we obtain 

IlSu! + 6w 2 {h -» (ft + 6w 3 ) (I 3 - I 2 » = 0 (C 4 ,a) 

“ 2^‘* ) 2 ~ 5 w^(h + ((2 + 6(03) (I3 - 1 ^)} ~ 0 (C 4 ,b) 


I36W3 = 0 (C 4 ,c) 

Thus we see immediately 60)3 - 0 , therefore W3 = 12 = constant. 
Then, equations a and b can be combined to yield 

6w x + (h + ft(I 3 - I 2 )}{h + n(I 3 - I 2 )> = 0 (C 5 ) 

I 1 I 2 


C-l 


d. 


m 1 


r 
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What we have obviously obtained here is the equation for a 
simple harmonic oscillator. Consider one special case first. If 
3 is "small", then ui (and W 2 * which has the same form of solution) 
oscillates with circular frequency 



'if a is not small, let us define an equivalent inertia 


I 3e = I 3 + h/a .(C7) 

Then eq. (C5) reads ; 

+ 6wia^(l3e ~ *2^ ^3e ” ^1^ = ^ (C8) 

It is clear that if either * 

l 3 e y I 1 a nd J 3e > *2 (C9,a) 

or 

I 3e < Ii and I 3e < I 2 (C9,b) 


that is, if l 3 e is a minimum or maximum inertia, then, our system 
is stable in the Liapunov sense. On the other hand, if l 3 e is an 
intermediate inertia, the system is unstable. These results are of 
a physically appealing nature since they are of the same form as 
the results for spinning rigid bodies. In fact, if in eq. (C?) 
h = 0 , l 3 e = I 3 , and those results are obtained directly. 

For non-rigid (i.e., damped) bodies, these requirements are 
modified somewhat. First, there is generally little damping on 
the rotor, and that which is present has a destabilizing effect 
(almost always) . Consider then a system with negligible damping 
on the rotor and a non-negligible amount on the main body. If 
the spin directions are such that the momenta of rotor and main 
body add (Q> 0 ) , it can be shown that it is necessary and sufficient 
for stability that l3e be the major principal axis (maximum) . Con- 
versely, if the momenta oppose each other (fl< 0 ) , it is necessary 
and sufficient that 

h Q « h + l 3 fl>0 (CIO) 

*-h_ <n<o (cii) 

13 

There is one additional region of stability for this system, and 
it exists only if I 3 is a major axis. The region defined by 

H < -h 

13-max (li,l 2 ) 


(C12) 





i 
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0 

3 


corresponds to l 3 e being a major principal axis while ft<0. As 
h —*0 the above results reduce to the simple major spin axis 
requirements. 

Definition (C7) may be combined with this observation into 
a simple rule, to wit 

h Q > in (CIO) 

where ho is the system momentum bias (equal to h + I 3 Q) , and I * 
max(Ii, 12)- This same conclusion is obtained in reference 6 
(page A- 82) . 

The result (CIO) has a direct application to the attitude 
acquisition phenomena noted under the section "Terminal Orientation" . 
Namely, since the magnets will drive the satellite at 2wo in the 
terminal orientation, and since we wish the terminal attitude to 
correspond to a stable condition, it is necessary that (CIO) hold 
for SI ~ 2uq* For a given orbit (which determines wo) and given 
set of inertias, this requirement then specifies a minimum h. The 
validity of this conclusion is supported by a number of simulation 
runs. In fact, it was observed that for good terminal attitude 
acquisition performance the left side of (CIO) should be at least 
three or four times the right side. 



» 

i 





it i > 
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APPENDIX D 


SYSTEM EQUILIBRIUM POINTS 


I 
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Consider the equilibrium states of the acquisition system 
described in this report under two conditions. (a) The magnetic 
field is fixed in trajectory ixis coordinates, that is 

d SPACE (g) = 0 (Dl) 

dt 

and (b) , the magnet controllers employ a linear law 

M * -KB (D2) 

• 

where B is the time rate of change of B as measured in the vehi'° 2 - 
As mentioned in the body of this report, under the section titled 
"Stability", this system has no equilibrium points unless some 
modification is made to the flip-flop magnet control law. No.. , 
let us consider the equilib ium points possible under the two con- 
ditions above. 

The behavior of our system may be completely characterized by 
the six state variables 012 , <*> 3 , Bi, B 2 , B 3 . Equilibrium points 

may be identified by setting the time derivatives of these variables 
to zero and solving the resulting six equations In six unknowns. 

That is, we require 
• 

tu — 0 (D3 ) 

• i 

B = B x a) = 0 eq. ( 8 ) — (D4 ) 

Starting with equation (1), we make the following substitutions: 
x * 0 torque free motion), H = 0 (constant speed wheel), and eq.(l') 

H = 5*5 + h. We then obtain 

1 * 0 ) = -u> x (i-w + h) = 0 (D5) 

The right hand side of (D5) equals zero according to (D3) . $n 
equilibrium point must then, by definition, satisfy equations (D4) 


and (D5) . (D4) simply requires that 

w = a(t)B (D 6 ) 

where a(t) may or may not be zero and may or may not be constant. 
Now, if 

h = {0, 0, h } = h (D7) 

That is, we have a pitch wheel only, then (D5) becomes, in scalar 
form 

W 2 U 3( I 3 ~ *2^ + = 0 (D 8 ,a) 

(D5, modified) 

wi« 3 (Il - I 3 ) - hwi = 0 (D 8 ,b) 

« 1 « 2( i 2 ** *l) *= 0 (D 8 ,c) 

D-l 




] 


£■ 


i 
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These four classes contain all the equilibrium points of this system. 
In passing I note (for the mathematicians among us) that there are 
an infinity of equilibrium points. Note also that at equilibrium 
points, § = o so M = 0 anc * so the f act that we set t “ 0 at the 

outset does not affect the validity of the results obtained here. 


J 

| From 

the above 

discussion we now obtain the following results j 


(I) 

w = (0 , 

1 

0, 0), B = any, a = 0 is an equilibrium ; 

? J 


point. 

since it clearly satisfies eq. (D4) & (D8). 


(II) 

» « (0, 

any, h/(I 2 “ 13 )}* B || w is an equilibrium ! 

j ] 

, 

point. 

i 

I 


(III) 

w = (any, 0, h/(Ii - X 3 ) > , B w is an equilibrium 

i 1 


point. 


V. 

i 

^ (IV) 

w - (0, 

0, any), B |[ u) is the last equilibrium point. 



There is, curiously, an alternate method of obtaining these 
same results. Starting with the equation of motion (1) , the control 
law (D2) , and the assumptions (Dl) and (D7) , we have 

ti = Iiwi + W 2 h (D9,a) 


T 2 = *2 a) 2 " w l^ 
t 3 = I 3 W 3 


where we have dropped the second order terms involving products of 


the u)’s. Using eq. (D 2 ) and 

(D 4 ) , we obtain 


= K(u>2B3 — U3B2) 


(D 10 , a) 

M2 — K(w3B3 — W3B3) 


(D 10 ,b) 

M3 = K(u)3B2 “ W2B1) 



0 

Remembering eq. ( 5 ) 


» 

T = M x B 

• 

( 5 ) 

and substituting t into eq. 

(D 9 ) , we obtain 


Iia»l + o)2^ 5=1 B 3 K(w3B^ - 

W1B3) - B2K(uiB2 “ W2 b x) 

(Dll, a) 

l20>2 ~ w^h = BiK(w^B2 ~ 

u 2 Bx) “ B 3 K(u> 2 B 3 ~ U3B2) 

(Dll,b) 

I3W3 = B 2 K(u> 2 B 3 - W 3 B 2 ) 

- BiK(u>3Bx - 

(Dll ,c) 


I 

1 
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Using the definition of equilibrium points 

w - o 

(here we are considering only <4 as a state variable, 
not B) 

we obtain three linear homogeneous equations in the three unknowns - 
t*>l, o>2, and <43. For a solution to exist, the determinant of the 
coefficients must equal zero. That is 


-(B 2 + b 3 ) 

(B^ - h/K) 

b 1 b 3 



(B3.B2 + h/K) 

-<B| + Bj) 

B? b 3 

= 0 

(D12) 

B 1 B 3 

b 2 b 3 

-»? + B 2> 




This equation reduces to 

-l4 (B l + B 2 } = 0 (D12 ' ) 

K 2 

or 

B = {0, 0, B 3 ) (D13) 

Substituting this back into eq. (Dll) we have 

+ u)2h = -B^Ko^ (D14 ,a) 

I 2 u 2 " “jh = -B 3 Ku 2 (D14 ,b) 

I3043 * 0 (Dl^,c) 

Since w = 0, we obtain the following equilibrium condition 

w = (0, 0, any) (D15) 

This result is not valid if U3 is large, since we neglected 

<43(13 - l 2 ) and <43(13 “ J l) -*- n comparison to h in writing eq. (D9). 

To the extent that this is valid, (D15) is correct. Note that (D15) 
is the same as condition (IV) identified previously, and that (I) 
is a special case of (IV). 

The only special thing about this alternate derivation is that 
we did. not assume g _ q at the outset. We were considering cases 
where ® = 0 but | and hence t not_necessarily equal to zero. Note 
that, in retrospect, we now have B = (o, 0, B3}, so 61 = B2 = Mi = 

M2 = 0. The result is then that t = 0 anyway. 

In the next section we consider the stability of the equilibrium 
points enumerated here. 
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APPENDIX E 


STABILITY OF EQUILIBRIUM POINTS 


t 




\ 
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In this appendix we investigate the stability of a particular 
equilibrium pcint, namely point IV of Appendix D. This stability 
investigation is carried out by linearizing the system equations 
about the equilibrium point and then forming the characteristic 
equation. A study of the roots of this equation then determines 
the stability of the system near the equilibrium point. This does 
not rule out the possibility of limit cycle like behavior but is 
nonetheless a useful technique. From the analysis that follows, 
certain conclusions may be drawn. These conclusions and their 
consequences are discussed at the end of this appendix. 

Consider perturbations about the equilibrium point 


a) = (0 , 0 , ft} 

B = { 0 , 0 , B} 

Then 

u> = '6(1)2, ft + 6(03} 

B — {6B1 , 6B2» B + 6B3} 

The following equations then apply 
B = B x ii) 

6Bq_ = ft5B2 “ B6 u>2 

6B2 = B6w^ — ft6B^ 

6B3 = 0 

where we have dropped second order terms in eq. 
a linear control law, i.e., 

M m -KB 

M X = K(B 5 u 2 - ft6B 2 ) 

M 2 = K(ft6Bi - BSwi) 

M3 = 0 

Then, since 

T = M x B 
We obtain 

Tl * KB(ft6Bx - B6wx) 

T 2 88 KB(ft6B2 “ B6 w2) 

T3 = 0 


(E3). 


(El) 


(E2) 


(8) 

(E3) 


Now, assuming 
♦ 

(E4 ) 


(E5) 


(5) 

(E6) 


E-l 
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Similarly, we may start with eq. ( 1 ), substitute, from (E 2 ) 
and drop second order terms in the result. Thus we obtain the 
linearized equations of motion about the equilibrium point 

I^6w^ = t ~ 6o> 2 (H + fl(l3 - I 2 )) (E7) 

I260J2 = ^2 + 6w^(h + n(I 3 - 1^) ) 

I36W3 = x 3 

These may be written in a way that makes them easier to interpret 
by using the following shorthand notation: 


I 3 e = h + I3 


(E8) 

ft 



I 3e ft “ h + I 3 JJ “ ho 


(E 9 ) 

h + ft (X3 - 12 ) - ftd 3 e “ 

12) = ai 

(E 10 ) 

h + ft (13 - I3.) = ft(! 3 e “ 

II) = a 2 

(Ell) 


Using this notation and substituting for the xi from (E6) we obtain 
60)! = -KB 2 6 o) ai6o) 2 + KBft6Bi + 06 B 2 

1 ' IT 

6u> 2 * 0261*)! - KB 2 6u ? + O6B1 + KBft 6 B 2 (E 12 ) 

TJ I 2 J 2 

6B 3 — +O60)! ~ B6 o)2 + O6B1 + Q6B2 
6B2 = B6oi! + 06 u >2 — 6 Bj_ + O6B2 

w 

The last two equations are from (E 3 ) . Note that the last equations 
in (E 3 ) and (E 7 ) immediately yield 

6B3 = 0 , 6B3 = const = 0 , B3 = B (E 13 ) 

60)3 = 0, 60)3 — const = 0, u> 3 = ft- 

Equations (E 12 ) constitute four first order ordinary linear dif- 
ferential equations in four unknowns. They are of the form 

X * A X 

where A is a square matrix with constant coefficients. 


(E 12 ,a) 
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Now, if A is negative (positive) definite it may be immediately 
shown that the system described by (El2,a) is assymptotically 
stable (unstable) by applying Liapunov's (Chetayev's) stability 
theorem to the function v ** x*x. Unfortunately, for our problem, 
as for many others, A is sign indefinite and the theorems noted 
are not so easily applied. 


Instead, we Laplace transform (El2,a) 

SX(S) - x 0 = AX (s ) (E14) 


Thus 


X(s) 

= (A - Is ) ~ 1 (“Xq ) 

and 


x(t) 

* L -1 { (A - IS)“ 1 (-X 0 ) } 

Note 

that z 

as shown by E15, all 


x(t) and X(s) are column vectors 
D(s) = Determinant (A - is) 


Now/ D(s) may also be written as 
D(s) = (s - si) (s - S2 )***(s 


(E15) 
(E16) 

the components of X(s) (remember, 
) have a denominator D(s) and 

(E17) 

- s n ) (E18) 


where the si are the roots of D(s) = 0 or Det (A - Is) «* 0. 
Alternately, the si are the eigenvalues of the matrix A. Obviously, 
for stability, it is necessary and sufficient that the real parts 
of the Si be less than zero. 


To determine the si we expand 


D(s) = Det (A - Is) = 0 


(E19) 


obtaining the (linearized system) characteristic equation in terms 
of the various system parameters. The matrix A from (E12) is 


A « 


H'-KB 2 /!! 
a2/ I 2 
0 
B 


-ai/Ii 

-kb 2 /i 2 

tB 


KBfl/12 

0 

0 

-n 


o 

KB0/I2 

ft 


and the expansion (E19) yields, after a little algebra 
Ill2S 4 +(Ii + I 2 )KB 2 s 3 + (Iil2« 2 + a x a2 + K 2 B 4 )s2 
+ KB 2 Q(ai + a2)s + aiq2fi 2 * 0 


(E20) 


(E21) 
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We return now to conditions (a3) and (a4). <{>i can always be 

made greater than rero for sufficiently large gain K. For K = 0 , 

(a3) may be written as 

2(11 + + (I 1 + 121 (I 3e - »3e - U) > 2I (E24) 

1112 36 

• 

Various cases of this can now be considered. If ft is small and 
less than zero (opposite in sense to h) , (E24) is easily satisfied. 

If ft is small and positive this relation is always satis f ad. For 
large values of ft (either sign) I3 e -c»l3. If I3 is a major or minor 
axis (E24) is always satisfied. If 13 is an intermediate axis, 

(E24) may fail (for K = 0) for sufficiently large ft (which is re- 
quired to make l3e*Ol3) • Before this happens the criterion (a5) 
will fail, so we need not concern ourselves with ( a 3 ) except perhaps 
for some intermediate values of ft, such that neither l3e >> Il» I2 nor 
I3 e — >13. If we make some simplifying assumptions- namely Ii* -• I2 = I 
and l3 e = nl, (E24) reduces to 

(n - l) 2 > (n - 2) (E25) 

which holds for all values of n, positive and negative. It would 
seem that except for a rather oddly configured satellite that 4>1>° , 
and (a3) presents no problems. In fact, unless the ratio between 
It and 12 is greater than 5:1, a3 must be satisfied. Additionally, 
for K>0, (j»i increases and (a3) presents even fewer problems. On the 
basis of this analysis we turn our attention to condition (a4) . 

$2 *» KB 2 ft(ai + 03) - KB 2 ft2a 1 o 2 (Ii + ^/^l > 0 ( a4 ^ 

<J> 2 = KB 2 ft 2 [(2I 3e -Ii-l2)-n 2 (Il+l2^ I 3e" I l ) (I 3e“l2)/<J)il >0 ( E26 > 

Now, from condition (a5) we see that the second term of (E26) is 
less than zero. Thus, if the first term is less than zero, this 
system is unstable. It is a requirement then, for stability, that 
l3e be the major principal axis of the satellite (ie, Ise^ax (Ii , I2) ) * 

Now, for our system, I3 is a minor axis, and for large values of 
8, l3e _ *' I 3/ so the system is unstable. For ft>0, in fact, I3e is 
always a minor axis (if I3 is). The only stable values for opera- 
tion are positive values r»«t "large”. In examining this problem 
we now turn to a numerical study. 

For a numerical evaluation, the following values are sub- 
stituted in conditions (a3) , (a4) and (a5) . 

I]_ “ 35 slug-ft 2 (47.5 Kg-m 2 ) 

12 * 50 slug-ft 2 (67.9 Kg-m 2 ) 

13 = 25 slug-ft 2 (33.9 Kg-m 2 ) 

B “0.5 Gauss 

h * 0.7 lbf-ft-sec (0.944 Nt-m-sec) 

K “ 10 4 pole-cm/(5 x 10" 4 Gauss/sec) *. 2xlO" 8 Nt-m/(Tesla 2 /sec) 



) 
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K is chosen by calculating the describing function for a flip-flop 
control law and a system oscillation amplitude of about uo* Note 
that K decreases as w grows and hence the system becomes less 
stable. Now, (a3) becomes 

♦l * 50 19 Q 2 - 97. 7« + 0.89 + K 2 B 4 > 0 (E27) 

(E27) has no real zeros even for the worst case (K = 0) , so we 
need not concern ourselves with it further. We next look at (a5) 
and determine the maximum value of f 2 for which is a major axis 
(recalling that this is required by 412) 

h + I3 > max(Ii, I2 ) (E28) 

. ft 

The result is ft max = 0,0278 rad/sec. 

Finally, we examine <$>2 (expression a4) : 

♦2 ■ - 53315ft 4 + 5185 ft 3 - (47.6 + 102.7 )ft 2 + 1.89ft (E29) 

KB 2 4,1 *1 <J>1 

We are interested in determining whether 4>2>0 for .0278>ft>0. This 
is required if we are to have any stable region of operation at all. 
This task is somewhat simplified by first examining <j>i. We see that 
4>1 is relatively constant for this range of ft. In fact, for K = 0, 

(<J>1 = 5015ft 2 - 97.7ft + 0.89) (E30) 

$l(C) = 0.89 

<j>l(.01) = 0.515 

♦ x (. 02) = 0.94 

4> 1 (.0278) = 2.05 

and 

^lmin = 0.418 at ft = 0.00973 

K = 0 is surely the worst possible case, since as K increases, so 
does <|>i and thus $2 will increase if we are in the range 0<ft<0.0278. 

Now, (E29) has only two real roots. One is obviously at 
ft = 0 , the other lies at about ft « 0.05. Between these roots 4*2 > 0. 
Thus we have ensured a stable region of operation for 

0.0278 > ft > 0.0 (E31) 

One may convince oneself of the truth of this result by re- 
turning to equation (E21) , substituting values for the coefficients 
and then solving the resulting equation numerically for the roots 
si. It is, of course, a necessary and sufficient condition for 
stability that all the si lie in the left, half plane. 


1 
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We require, for stability, that all the solutions of (E21) lie in 
the left half-plane. We may attack (E21) either numerically or 
analytically. Consider the analytical approach first. 

If we apply Routh's Criterion to (E21) we obtain, as the 
necessary and sufficient conditions for stability, that the 
following quantities must le greater than zero if ai is 


Ill2 > 0 (al) 

(II + I 2 )KB 2 > 0 (a2) 

4»1 = Iil 2 ^ 2 + aict 2 + K 2 B^ - + a2) > 0 (a3) 

II + 1 2 

<J >2 = KB 2 f2(ai + a2) - KB 2 n 2 aia 2 (Ii + l 2 )/<|>]. > 0 (a4) 

a l a 2^ 2 > 0 (a5) 


Requirement (al) is trivial, and so is (a2) as long as K>0. (a5) 

is satisfied as long as the body is dynamically stable (see appen- 
dix C) without the control system. We must address ourselves to 
(a3) and (a4). We first note, en passant, that if we set K = 0 
in (E21) we obtain the characteristic equation for the system with- 
out control torques. This situation was examined in appendix C. 
Equation (E21) becomes 

s 4 + as 2 b = 0 (E22) 

where 

b = aic^ft 2 
~Ill2 


and 


a = n 2 + ai «2 + K 2 B 4 

I]d2 

It is trivial to show that the requirements for stability of (E22) 
are that a>0, b>0. Obviously b>0 is the stronger condition; we 
obtain the requirement 

a x a 2 > 0 (E23) 

and this is always satisfied if the system is stable without the 
control system. 


* 


-N 


# 


! 
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Substituting into (E21) we obtain 

3225s 4 + 57.7s 3 + (3687a 2 - 45a + 1.14)s 2 + (-23. 8Q 2 + 0.944.°.)s 

+ (462a 4 - 45a 3 + 0.89Q 2 ) = 0 (E32) 

Numerical solutions for the roots of (E32) confirm the results 
(E31) obtained thru the Routhian analysis. 

All of the results obtained in this appendix# up to this point, 
are only valid in a rigorous sense, for motion in the "neighborhood" 
of the equilibrium point. Otherwise the linearizations made in 
obtaining the system equations (El 2) are not allowable. Nonetheless, 
considerable digital simulation experience indicates that the res Its 
obtained in this section are in fact applicable to the real (noi 
linear) system. Let us consider then the significance of these re- 
sults. First of all, they indicate that the terminal attitude of the 
satellite must be such that the momentum bias points along the orbit 
normal, rather than opposite it. We remember that in the terminal 
condition the satellite has an angular velocity 2oj 0 e n . If the satel- 
lite were to acquire "backwards", then the angular velocity 2uoe n 
would correspond to a<0, and we have seen that this equilibrium point 
is an unstable one. We have seen that there are no stable equilibrium 
points with very high angular velocities; the greatest being a = 

0.0278 rad/sec. This means that, with the exception of possible limit 
cycle difficulties, this system should be well behaved and the control 
theory developed in the main body of this report provides for the 
positive despin of the vehicle. The author also wishes to note that 
he has observed no limit cycle behavior during extensive digital simu- 
lation of this system. 

The reader is probably asking, at this point, "How about the 
equilibrium points II and III enumerated in Appendix D" . The answer 
is that the author has not yet had time to perform a rigorous stabil- 
ity ana 1 ysis for these points. However, on the basis of digital simu- 
lation rt appears that they are not stable points unless ui = w2 ~ 0, 
which is of course the case treated here. 

Finally, we note that even the stable equilibrium states deter- 
mined here present no real despin problems. The investigations here 
dealt, for simplicity, with a fixed external magnetic field. In fact, 
even in the worst case (an Equatorial orbit with the pitch wheel along 
the orbit normal) the Earth's field appears to cone about with a half 
angle of 11.7°. If the satellite is to remain in the equilibrium 
state, it too must cone about the orbit normal at frequency u 0 * The 
control torques maintaining the equilibrium state (i.e. damping out 
disturbances - which is how the coning of the B field appears to the 
system) are essentially dissipative in nature, since they oppose 
sensed velocities. The final result is then that even the small 
residual spin f < .0278 rad/sec along the momentum bias will be elim- 
inated due to this previously unconsidered effect. 




